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Tuneable  Interferometric  Bandpass  and  Bandstop 
Filters  for  Terahertz  Applications 

Otto  Schwelb  and  Istvan  Frigyes 


Abstract  - Numerically  simulated  performance  of  a new 
type  of  tuneable  bandstop  and  bandpass  filter  is  presented. 
The  filters  use  2x2  couplers  loaded  by  discontinuity-assisted 
ring  resonators  (DARR).  The  bandstop  filter  consists  of  one 
or  several  cascaded  DARRs.  Two  bandstop  filters  separated 
by  a cavity  make  up  the  bandpass  filter.  Tuning  is  achieved 
by  varying  the  geometrical  or  material  parameters  of  the 
ring  resonators  or  the  cavity.  Alternatively,  when  the  device 
is  used  as  a sensor,  environmental  conditions  are  monitored 
by  their  detuning  effect.  Selectivity  and  tuning  range  are 
controlled  by  the  coupling  coefficient  of  the  couplers,  by  the 
circumference  of  the  rings  and  by  the  reflection  coefficient 
of  the  discontinuities.  The  effects  of  coupler  and  waveguide 
loss  have  also  been  investigated. 

I.  Introduction 

Quasi-optic  tuneable  bandstop  and  bandpass  filters  with 
bandwidths  from  the  KHz  to  GHz  range  are  key 
components  in  measurement,  telecommunication  and 
sensing.  Analytical  investigations  show  that  a 2x2 
coupler  (lumped  element  or  distributed  parameter  type) 
loaded  by  a ring  resonator  that  incorporates  a controlled 
amount  of  small  reflective  discontinuity  can  be  a very 
effective  bandstop  filter  subject  to  the  strength  of 
coupling.  Using  this  circuit  as  a basic  building  block 
(BBB)  we  investigated  several  bandstop  and  bandpass 
filter  configurations  with  remarkably  selective  properties. 
The  bandstop  filter,  effectively  a highly  reflective 
composite  mirror,  is  fabricated  from  one  or  several 
cascaded  BBBs.  The  bandpass  filter  consists  of  two  such 
mirrors,  separated  by  a cavity,  much  like  a Fabry-Perot 
resonator.  The  discontinuity  built  into  the  ring  resonator 
can  be  a mismatch  or  a short  length  of  waveguide  with 
material  or  geometric  properties  slightly  different  from 
that  of  the  ring. 

The  results  have  been  obtained  from  a computer  program 
using  a lumped  element  model  for  the  couplers;  neither 
the  fabrication  technology  of  the  resonant  ring  nor  that  of 
the  discontinuity  is  addressed  since  they  strongly  depend 
on  the  band  of  operation.  Some  of  our  simulations  have 
been  carried  out  at  3THz,  others  at  194THz. 

Tuning  of  the  bandstop  filter  is  achieved  by  modulating 
the  parameters  of  the  ring  or  that  of  the  discontinuity. 
Tuning  of  the  bandpass  filter  requires  variation  of  the 
phase  delay  in  the  cavity  between  the  mirrors.  To  obtain 
such  a phase  delay  variation  either  the  length  or  the 
refractive  index  of  the  cavity  must  be  adjusted.  Among 
the  means  to  produce  this  adjustment  are 
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micromechanical,  piezoelectric,  electrooptic  and  thermal. 
Conversely,  environmental  conditions  can  be  monitored 
by  their  detuning  effect.  The  technologies  of  tuning  and 
modulation  have  not  been  here  considered. 

11.  Discontinuity-assisted  ring  resonator 

The  BBB  of  the  filter  is  a 2x2  coupler  loaded  by  a DARR 
as  shown  in  Fig.l,  where  K is  the  power  coupling 
coefficient  of  the  coupler,  y its  fractional  power  loss  and  r 
is  the  reflection  coefficient  of  the  discontinuity. 


Fig.l.  The  basic  building  block  of  the  filters:  a 2x2 
coupler  loaded  by  a discontinuity-assisted  ring  resonator. 

It  is  apparent  that  without  the  reflective  discontinuity  the 
circuit  has  an  all-pass  characteristic  [1],  i.e.,  the 
magnitude  of  the  transmitted  intensity  in  a lossless  device 
is  unity,  and  therefore  without  the  reflective  discontinuity 
the  lossless  circuit  cannot  be  used  as  a frequency  filter. 
(A  lossy  circuit  without  a discontinuity  will  still  not 
reflect  power  back  to  the  input,  but  its  transmission  will 
exhibit  a notch  near  resonance.)  Indeed,  by  exploiting  the 
frequency  dependency  of  its  transmission  phase,  this 
circuit  can  be  used  as  a group  delay  compensator  [2]. 
While  its  reflectivity  must  be  carefully  controlled,  the 
nature  of  the  discontinuity  embedded  in  the  resonant  ring 
is  immaterial.  It  can  consist  of  a small  mismatch,  a single 
step-index  dielectric  layer,  a short  shallow  grating,  etc. 
We  assume  a symmetric  discontinuity,  which  is  most 
conveniently  represented  by  the  scattering  matrix  [3] 


Jt  r J 


where  t = yj\-r^  is  the  transmission  coefficient  of  the 
discontinuity.  In  practice  the  reflection  coefficient  is 
small,  in  the  order  of  10'''  to  10''  therefore,  to  obtain  a 
high  Q resonance,  the  power  coupling  coefficient  K must 
be  small  as  well,  usually  between  40dB  and  lOdB.  Note 
also  that  it  makes  no  difference  for  the  operation  of  the 
DARR  where  inside  the  ring  the  discontinuity  is  located. 
Using  a conventional  model  for  the  coupler  [4]  the  two 
distinct  elements  of  the  scattering  matrix  of  the  BBB  are 
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-a  rK  exp(  jkL) 

5,  = ; r-= (2) 

exp{j2kL)  -a  K)-  j'2at\\  - K exp{jkL) 
and 


[exp(  /2AL)  - xjl-  K - jat(2  - K)  exp(jkL) 

~ ^ 

exp{j2kL)  -a  (I  - K)~  j2at\\-  K e.xp{jkL) 

where  k is  the  wavenumber,  L is  the  physical  length  of 
the  ring  and  a = is  the  amplitude  loss  coefficient 

of  the  coupler  (;ris  the  fractional  power  loss). 

Equations  (2)  and  (3)  reveal  that  the  BBB  is  a bandstop 
filter  whose  behaviour  is  subject  to  the  value  chosen  for 
K.  Fig.2  depicts  the  wave  intensity  transmitted  through 
the  BBB  as  a function  of  detuning  from  the  design 
frequency  and  as  a function  of  K.  (To  convert  from 
detuning  to  frequency  multiply  by  c the  free  space 
velocity  of  light  in  mm/s.)  The  figure  shows  that  at 
design  frequency  there  is  a critical  value  for  K,  denoted 
by  K„  where  the  transfer  intensity  is  zero.  For  K < the 
resonator  is  undercoupled,  showing  two  minima, 
increasingly  separated  as  the  value  of  K decreases.  For  K 
> Kc  the  resonator  is  overcoupled  showing  one  minimum 
which  becomes  increasingly  shallow  as  the  value  of  K 
rises.  The  relationship  between  K^.  and  r is: 

2r 

K^= s2/-  (4) 

1 + r 

Important  to  notice  is  the  second  derivative  of  Iz\  with 
respect  to  frequency  at  centre  wavelength.  For  an 
undercoupled  DARR  this  is  negative,  whereas  for  an 
overcoupled  DARR  it  is  positive. 


Fig.2.  Intensity  transfer  through  the  BBB  as  a function  of 
detuning  and  power  coupling  coefficient. 

111.  BANDSTOP  FILTER  CONSTRUCTION 

Although  the  BBB  offers  large  rejection  levels,  its  skirt 
selectivity  is  not  sufficiently  steep.  This  can  be  improved 
upon  by  cascading  several  BBBs,  separated  by  lengths  of 
waveguides.  To  reinforce  the  reflection  from  each  unit, 
the  phase  delay  between  consecutive  BBBs  must  be  such 
as  to  cause  constructive  interference  at  the  entrance  port. 


much  like  in  the  case  of  reflected  wavelets  from  each 
element  of  a Bragg  grating  in  the  stopband  region  [5]. 
Such  constructive  interference  is  achieved  when  the 
length  of  the  waveguide  separating  the  BBBs  is  an 

integral  multiple  of  = A^2n^,  where  To  is  the  design 
wavelength  and  the  effective  index  of  the  guide.  Fig. 3 
shows  the  transfer  characteristics  of  N cascaded, 
undercoupled,  identical  BBBs.  The  bandwidth  of  the 
stopband  of  this  composite  mirror  is  proportional  to  r and 
inversely  proportional  to  L^jj-  the  effective  length  of  the 
resonator  (lengthxeffective  index): 

cr 

Af  = (5) 

This  is  in  accordance  with  the  fact  that  the  group  delay 
through  the  composite  mirror  at  To  is  proportional  to  L^jj 
and  inversely  proportional  to  r.  The  performance 
characteristics  of  the  composite  mirror  in  the  stopband 
region,  including  the  group  delay,  are  similar  to  those  of  a 
long  Bragg  grating  [6].  By  cascading  undercoupled  and 
overcoupled  building  blocks  square  well  (flat  bottom) 
stopband  characteristics  can  be  obtained. 


Fig.3.  Intensity  transfer  as  a function  of  detuning  through 
jV  cascaded,  undercoupled,  identical  BBBs. 

IV.  Fabry-perot  type  bandpass  filter 

Since  the  transmission  characteristic  of  a cascaded  BBB 
is  analogous  to  a composite  mirror,  one  can  build  a 
Fabry-Perot  type  bandpass  filter  by  inserting  a waveguide 
between  two  identical  chains  of  BBBs  as  shown  in  Fig.4. 


Fig.4.  Fabry-Perot  type  bandpass  filter  consisting  of  two 
composite  mirrors  separated  by  a cavity. 

If  the  length  of  the  cavity  (or  defect,  as  it  is  sometimes 
called)  is 
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A 

=(«.»  + integer  (6) 

then  the  passband  peak  will  fall  in  the  middle  of  the 
stopband  of  the  composite  mirrors,  i.e.,  on  the  design 
frequency.  The  maximum  continuous  tuning  range  within 
the  mirror  stopband,  obtained  by  varying  the  cavity 

length  is  approximately  AL^av  = O.Sf^,  where  AL,.a^,  is  the 

change  in  the  physical  length.  Larger  values  of  AL,.„v  to 
obtain  the  same  tuning  range  can  be  realised  by  jumping 
from  one  n^av  to  another,  as  implied  by  (6).  When  the 
cavity  is  tuned  by  refractive  index  variation  the  maximum 
tuning  range  is  An,.  = «t./(2«,.„„+l).  In  either  scheme  the 
total  phase  delay  while  tuning  the  cavity  from  end  of  its 
range  to  the  other  is  approximately  A^=nl2.  Fig.5  shows 
the  intensity  transfer  characteristics  of  the  filter,  centred 
at  3THz,  using  composite  mirrors  built  with  three  unit 
cells  having  coupling  coefficients  Ai,=Ai3=0.00103 
(29.9dB)  and  /f2=0.00255  (25.94dB),  the  reflection 
coefficient  was  r = 0.001.  The  solid  lines  are  the 
bandpass  curves  corresponding  to  cavity  lengths  L,-a„  = 

100.58^,1,  100.504  and  100.424,  respectively,  the  dash- 
dotted  line  marks  the  flat-floor  bandstop  of  the  mirrors. 
The  passband  peaks  rise  from  this  floor;  the  lower  the 
floor  level  the  narrower  are  the  passbands.  A flat  floor 
insures  a constant  bandwidth  across  the  tuning  range. 
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Fig.5.  Transmission  intensity  as  a function  of  detuning 
for  three  cavity  lengths.  The  dash-dotted  curve  refers  to 
the  bandstop  filters  that  comprise  the  mirrors. 

V.  The  EFFECT  OF  LOSS 

The  finite  bandwidth  and  the  corresponding  finite  Q of 
the  bandpass  filter  fabricated  with  lossless  components  is 
an  indication  that  the  composite  mirrors  are  partially 
transmitting.  The  smaller  this  transmission  is,  the  larger  is 
the  Q: 

Q = 27t-^  (7) 

where  /21  is  the  numerical  value  of  the  floor  level 
transmission.  Losses  reduce  the  peak  transmission  and 
broaden  the  resonance  curve. 


In  interferometric  circuits  one  encounters  coupler  loss 
and  ring  resonator  waveguide  loss.  These  have 
demonstrably  equivalent  effect  on  the  transmission  and 
group  delay  characteristics  of  the  filter.  In  comparison, 
waveguide  losses  in  the  connecting  guides  and  in  the 
cavity  have  insignificant  effect  on  transmission.  This  is 
because  a resonant  ring  'magnifies'  its  losses  (as  well  as 
its  transmission  group  delay)  proportional  to  MK,  usually 
a large  number.  Fig.6  is  a surface  plot  of  the  bandpass 
peak  in  the  vicinity  of  band  centre  as  a function  of 
fractional  coupler  loss.  The  plot  shows  the  rapid 
deterioration  of  the  passband  characteristic  with 
increasing  y.  The  parameters  used  to  obtain  Fig.6  are  the 
same  as  those  for  Fig.5.  An  identical  plot  is  obtained 
when  instead  of  increasing  y from  0 to  y„ax,  the 
waveguide  attenuation  coefficient  a is  increased  from  0 
to  a^ax  = ymaJ'2.1,  whoro  a is  in  Np/mm  and  the  length  of 
the  ring  is  in  mm.  Respecting  the  above  law  of 
equivalence  between  attenuation  coefficient  and 
fractional  power  loss,  the  two  loss  mechanisms  are 
additive. 


Fig.6.  Surface  plot  of  the  intensity  transfer  as  a function 
of  detuning  and  fractional  coupler  loss. 

Transmission  levels  of  lossy  filters  also  depend  on  A'/Ai,. 
and  the  number  of  BBBs  in  the  composite  mirrors.  This  is 
illustrated  in  Fig.7,  showing  the  transmitted  intensity  of 
bandpass  filters  at  3THz  design  frequency  as  a function 
of  the  relative  coupling  coefficient.  The  solid  (dash- 
dotted)  lines  refer  to  y=  1x10'^  (;i'=  IxlO'A,  A represents 
the  number  of  identical  BBBs  in  each  mirror  and  or  = 0. 
Observe  how  the  transmission  level  for  a given  setting  of 
y drops  as  K approaches  K,.  and  how  sensitive  this 
transmission  is  to  N.  The  acute  sensitivity  of 
interferometric  circuits  based  on  resonant  rings  to  loss 
might  require  inclusion  of  active  elements  that  can 
provide  the  small  amount  of  gain  necessary  to  offset  the 
losses. 

VI.  Design  PROCEDURE 

The  fundamental  parameters  and  technological 
constraints  given  when  designing  a tuneable  ring 
resonator  based  terahertz  filter  are:  the  centre  wavelength, 
the  tuning  range,  the  stopband  or  passband  width,  the 
effective  refractive  index  of  the  guides  and  the  minimum. 
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or  desirable  length  (circumference)  of  the  rings.  The 
minimum  free  spectral  range  (FSR)  might  also  be 
specified,  but  its  value  is  subject  to  the  minimum  ring 
radius. 


Fig.7.  Transmission  intensity  vs.  relative  coupling 
coefficient  of  lossy  bandpass  filters  using  N identical 
BBBs  in  each  mirror. 


via  index  modulation,  increasing  the  value  of  n^av,  as 
indicated  above,  can  reduce  the  modulation  depth 
necessary  to  tune  through  the  stopband.  The  transmission 
and  reflection  characteristics  of  the  bilaterally  symmetric 
Fabry-Perot  type  bandpass  filter  are  obtained  from  its 
scattering  parameters; 


c _ c 

^/i/'ll  .\/ll 


1+- 


(12) 


and 


- 


(13) 


where  S^u  and  are  the  scattering  parameters  of  the 
composite  mirrors  on  either  side  of  the  cavity.  A 
computer  code  to  follow  through  these  steps  has  been 
developed. 


Vll.  Conclusions 


From  (5)  one  can  determine  the  reflection  coefficient  of 
the  discontinuity  r.  When  the  discontinuity  is  a single 
defect  layer,  its  refractive  index  % and  length  4 are 
determined  from 


r = l-exp(-y2^ /J 
^ I- p'  exp(-j2kJJ 


(8) 


where  ka  is  the  wavenumber  in  the  defect  layer  and 


n + n ,, 

a 


(9) 


Simulated  characteristics  of  a novel,  tuneable  bandstop 
and  bandpass  filter  have  been  described.  The  bandstop 
filter,  fabricated  with  one  or  several  cascaded  2x2 
couplers  loaded  by  discontinuity-assisted  ring  resonators 
(DARR)  is  regarded  as  a composite  mirror.  The  bandpass 
filter  consists  of  two  such  mirrors  in  a Fabry-Perot 
configuration.  The  bandwidth  is  adjustable  and  can  be 
made  arbitrarily  small.  The  effect  of  the  parameters 
controlling  the  tuning  range  and  the  passband  width,  as 
well  as  the  effect  of  waveguide  and  coupler  loss  has  been 
investigated.  A detailed  design  procedure  has  been 
described.  The  technological  aspects  of  device  fabrication 
and  tuning  have  not  been  addressed. 
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assuming  »a  > n^ff.  The  maximum  value  of  r provided  by 
(8)  is 


\ + p' 


(10) 


and  the  relative  index  and  length  of  the  defect  layer  are, 
respectively 


n 1 + r , . I- r 

_2_=  and  / =(2A  + l)^^  =2;^  (11) 

n,  Vl-r  4n,\l  + r 

With  the  stopband  width  and  the  design  wavelength  we 
first  determine  the  Q,  then  substituting  the  value  of  r and 
Q into  (7)  the  floor  level  transmission  intensity  of  the 
square  well  bandstop  characteristics  is  computed.  Next 
the  coupling  coefficients  are  adjusted  to  meet  this 
requirement.  The  number  of  DARR  loaded  couplers  is 
determined  by  the  required  filter  skirt  selectivity  (and  loss 
considerations),  while  the  length  of  the  guide  separating 
the  couplers  within  the  composite  mirror  must  be  an 

integral  multiple  of  f The  length  of  the  cavity  at  centre 
frequency  is  determined  by  (6).  When  the  cavity  is  tuned 
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